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We show how methods of continuum perturbation theory can be used to simplify perturbative lattice calcula¬ 
tions. We use the technique of asymptotic expansions to expand lattice loop integrals around the continuum limit. 
After the expansion, all nontrivial dependence on momenta and masses is encoded in continuum loop integrals 
and the only genuine lattice integrals left are tadpole integrals. Using integration-by-parts relations all of these 
can be expressed in terms of a small number of master integrals. Four master integrals are needed for bosonic 
one loop integrals, sixteen in QCD with Wilson or staggered fermions. 


1. Introduction 

In many cases, the results of numerical simula¬ 
tions of QCD on the lattice need to be matched 
by equally precise perturbative calculations in lat¬ 
tice regularization to become phenomenologically 
relevant. In the past, the need (or the motiva¬ 
tion) for precise matching calculations has been 
limited by the large systematic errors induced by 
the quenched approximation. By now, simula¬ 
tions with fairly light dynamical quarks are fea¬ 
sible and precise calculations in lattice perturba¬ 
tion theory are becoming increasingly important. 
The dominant uncertainty in the first determi¬ 
nations of the strong coupling constant from lat¬ 
tice simulations with three dynamical quarks, for 
example, are unknown higher order perturbative 
corrections to the relation between the coupling 
measured on the lattice and the MS coupling [1,2], 

The introduction of a space-time lattice to reg¬ 
ulate the theory considerably complicates pertur¬ 
bative calculations. The propagators and vertices 
become functions of sines and cosines of the mo¬ 
menta and not a single loop integral in four di¬ 
mensions can be evaluated analytically. In the 
standard approach to lattice perturbation theory 
the relevant diagrams are therefore evaluated nu¬ 
merically, which has several drawbacks: (a) the 
amount of numerical computations necessary for 
realistic calculations is huge; (b) cancellations be¬ 


tween individual diagrams can render numerical 
results unstable; (c) the continuum limit, i.e. the 
limit in which the inverse lattice spacing becomes 
much larger than external momenta and masses, 
has to be taken numerically as well. There are 
a number of techniques to both reduce the num¬ 
ber and increase the precision of the numerical 
integrations involved [3-5], but we feel that the 
computational tools for perturbative calculations 
on the lattice are not as highly developed as the 
methods used in the continuum. 

In a recent paper [6], we have proposed to 
use continuum methods to simplify lattice cal¬ 
culations. In particular, two tools originally de¬ 
veloped for the evaluation of continuum loop 
integrals turn out to be useful: the technique 
of asymptotic expansions [7] and the use of 
integration-by-parts relations [8]. The method of 
asymptotic expansions can be used to expand lat¬ 
tice loop integrals around the continuum limit. 
This expansion splits the loop integrals into two 
parts: the “hard” part is a sum of lattice tad¬ 
pole integrals, while the “soft” part consists of 
ordinary continuum integrals. In a second step, 
we use integration-by-parts relations to express 
all lattice tadpole integrals through a few mas¬ 
ter integrals. In a theory involving only bosonic 
fields, there are four such master integrals at one 
loop. In QCD with Wilson or staggered fermions, 
the number of master integrals increases to six¬ 
teen. The master integrals can be chosen to be 
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convergent and are evaluated numerically. 

We begin by quickly reviewing perturbation 
theory in lattice regularization and some of the 
complications that arise with this particular reg¬ 
ulator. After illustrating the above method with 
a simple example, we discuss its application to 
the case of staggered fermions. 


2. Lattice perturbation theory 


When calculable at all, phenomenological re¬ 
sults in QCD are often obtained in a factorized 
form: a physical quantity is given as a product 
of a perturbative short-distance part times a low- 
energy contribution. The splitting into the two 
parts is regularization dependent. If one man¬ 
ages to evaluate the low energy part with a lattice 
simulation, one also needs the high energy part in 
the same regularization. This requires a pertur¬ 
bative calculation in lattice regularization. Such 
calculations are, however, tedious. First of all, 
the propagators are much more complicated than 
in the continuum. The simplest discretization of 
the bosonic propagator in four dimensions is 


G B (k) = 


(k 2 + m 2 ) 


(1) 


where 


k 2 = k„ and k tl = 2 si 
m= i 
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of gluons. The four-gluon vertex arising from the 
simplest Wilson loop, for example, involves one 
hundred terms with up to four powers of sines 
and cosines of the in- and outgoing momenta (see 
e.g. [9]). 


3. Asymptotic expansion around the con¬ 
tinuum limit 


We now discuss the expansion of the massive 
tadpole integral 


G(m) 


r d 4 k 1 

J-TT ( 2 tt ) 4 ( k 2 + m 2 ) 


(4) 


around the continuum limit, corresponding to the 
expansion around m = m p h ys a = 0. Altough this 
is the simplest possible example, it illustrates all 
the important features of our approach. 

We begin by mapping the integration region 
in Eq.(3) to an infinite volume and defining new 
integration variables 77^, 


Vu = tan(fc /J /2). 


(5) 


In terms of the new variables, the loop integra¬ 
tions in Eq.(4) range from —00 to +00: 


G(m) 


1 f TT dijj to 2 

471-4 J f 1 t 1 + ^ I- 4 

-OO 4-1 


Db(v) 

( 6 ) 


where 


Loop integrals involve products of these propaga¬ 
tors, integrated over the Brillouin zone; a typical 
loop integral has the form 


7T 



— 7T 


d 4 k 1 _ 1 

(2tt) 4 ( k 2 + m 2 ) ((p + k ) 2 + m 2 ) 


(3) 


where p is the external momentum. Note that we 
have expressed all quantities in units of the lattice 
spacing a: = fcp hys a, m = m phys a, etc. 

In lattice gauge theories not only the propaga¬ 
tors, but also the vertices take a rather compli¬ 
cated form. The gluonic action is given in terms 
of Wilson loops, whose expansion in terms of the 
gauge field yields interactions among any number 


Db{v) 


yy 'li 

i—l (1 + Vi) 


(7) 


After the variable transformation, the form of the 
integrand becomes more similar to a continuum 
integral. In addition to the change of variables, 
we have also introduced an intermediate regulator 
by raising the propagator to the power 1 + 5 in 
Eq. (6). 

The expansion of G(m) is obtained as a sum of 
two parts: 


G(m) = G so f t (to) + G hai . d (m). (8) 


The soft and hard contributions are calculated by 
applying the following procedure to the integrand 
in Eq.(6): 



3 


• Soft: assume that all the components of the 
loop momentum 77 are small, 77 ,; ~ m <C 1. 
Perform the Taylor expansion of the inte¬ 
grand in Eq.( 6 ) in the small quantities 77 , 
and m. The expansion coefficients in this 
region are standard continuum one loop in¬ 
tegrals, regularized analytically. Note that 
no restriction on the integration region is 
introduced. 


• Hard: assume that all the components of 
the loop momentum are large, 77 ,; ~ 1 777 

and Taylor expand the integrand in m. The 
expansion coefficients are massless lattice 
tadpole integrals. 

As in dimensional regularization, scaleless inte¬ 
grals are set to zero. Higher orders of the soft part 
contain more and more ultraviolet divergent inte¬ 
grals. The ultraviolet divergences manifest them¬ 
selves as poles in <5: 
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(9) 


In the final result, these poles will cancel against 
infrared divergences appearing in the expansion 
of the hard part, whose evaluation we discuss in 
the next section. 


4. Integration-by-parts relations between 
lattice integrals 

While the soft part is given in terms of contin¬ 
uum integrals, the hard part involves the func¬ 
tions 


00 4 

H (Wih n ) = / n (1 _^2 yl lDB(v)r n - d ■ 

-00 * =1 

( 10 ) 

The functions H({ai}; n), with aq and n integers, 
provide the full set of genuine lattice integrals 
needed for any one-loop calculation. The inte¬ 
grals occuring in the hard part are process in¬ 
dependent and can therefore be calculated once 


and for all for a given lattice action (the hard in¬ 
tegrals occuring in the case of staggered fermions 
are given in section 5). 

The hard part of Eq. ( 6 ) is obtained by expand¬ 
ing the integrand in m 

2 

4t r 4 G(m) hard = H( 1,1)-^- (1 + 6) H{ 1, 2) 

+ (1 + <>) ( 2 + S) H( 1,3) + ... , 

with 1 = {1,1,1,1}. The hard part of integrals 
with external momenta or tensor structure also 
involves functions H with index values a,; 7 ^ 1. 

We now use integration-by-parts identities to 
fully exploit the algebraic relations between the 
various integrals ff({a,;}, n). These relations are 
derived using the fact that the integral of a total 
derivative vanishes in analytic regularization: 



1 

{1+Vi) ai 


lD B (vT n - s }, 


( 11 ) 


for each value of 77 = 1,2, 3,4. Computing the 
derivative in Eq.(ll) and rewriting the resulting 
expression in terms of the functions H({ai},n), 
we obtain an algebraic relation between integrals 
with different values of { 07 } and n. Another equa¬ 
tion can be obtained by partial fractioning, i.e. by 
using the linear dependence of the five “propaga¬ 
tors” in the function H({a,i},n). The complete 
set of algebraic relations is therefore 

0 = {n'+EK + -l)k ( 12 ) 

^ i= 1 ' 

0 = |l + 2oj (af — 1) 

+2 (?r + 6) n + al(af - l)ji4. 

The conventions are such that the operator af 
increases (decreases) the index a,: by one. 

Similar integration-by-parts relations for lat¬ 
tice integrals were first studied in [3,4], where 
it was shown that the entire class of integrals 
H({ak}',n) can be reduced to d master integrals 
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in d dimensions. Here, we neither attempt to 
solve these equations explicitly nor to rewrite 
them in such a form that the reduction of a given 
index is manifest. Instead, we adopt a brute force 
strategy and use computer algebra to explicitly 
solve the equations for a given range of indices. 
An efficient algorithm for solving such recurrence 
relations has been described in [10]. First, a cri¬ 
terion which selects a simpler integral out of any 
two integrals is chosen. Typically, integrals with 
lower values of the indices are considered to be 
simpler. The above equations are then solved 
for a very limited range of indices, using Gauss’s 
elimination method. The calculation is repeated 
after supplementing the chosen set of equations 
with a few relations involving higher index values. 
By iterating this procedure, the equations (12) 
can be solved for the entire index range needed in 
a given calculation. The advantage of this brute 
force method is that it immediately generalizes to 
integrals involving more complicated propagators 
(e.g. those of Wilson fermions) or to higher loops. 

It is possible and advantageous to choose very 
convergent master integrals, e.g. the integrals 



with n = 0,1,2,3. The expansion of these in¬ 
tegrals in 5 is then obtained by expanding the 
integrand and it is trivial to numerically evaluate 
them to arbitrary precision. The tadpole integral 
with two powers of the boson propagator, for ex¬ 
ample, is 
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The result for the hard part of the massive tad¬ 
pole in Eq. ( 6 ) has the form 

G h *"V> = fa +™ 2 (-fa + +1-4)) 

+ ”' 4 ( ! ”-T2 bQ- hl4 )) +c,(m6> ' 

Note that the hard part is an analytic function 
of m. All nonanalytic dependence on masses and 
momenta arises in the soft part. As mentioned 
earlier, the 1/5-poles cancel between the hard and 
the soft part. The constants b\, 62 , 63 are ob¬ 
tained by numerically evaluating the master inte¬ 
grals. It turns out that these three constants are 
sufficient to express the hard part of any bosonic 
one-loop integral to any order in the expansion 
around the continuum limit. Their numerical val¬ 
ues are given in [ 6 ]. 


5. Staggered fermions 


The strategy outlined in the previous two sec¬ 
tions did not rely on the specific form of the prop¬ 
agator. In [ 6 ], we have applied the technique to 
HQET and to QCD with Wilson fermions and cal¬ 
culated a number of QCD one-loop self-energies. 
The recursion relations between the hard inte¬ 
grals are more complicated than in the bosonic 
case and the number of master integrals is larger. 
There are seven master integrals for HQET and 
sixteen for QCD with Wilson fermions. 

We now show how to use the method to calcu¬ 
late loop integrals for staggered fermions, which 
are a common choice to put fermions on the lat¬ 
tice. The staggered fermion action is obtained 
by reducing the number of components of the 
fermion field after spin diagonalizing the naive 
lattice fermion action. In perturbative calcu¬ 
lations it is convenient to work with the naive 
fermion action and perform the staggering only 
at the end of the calculation. The fermion dou¬ 
bling inherent in the naive discretization of the 
fermion action manifests itself in the appearance 
of multiple soft regions. The sixteen zeros of the 
propagator denominator of naive fermions, 


Df = T 7 sin 2 ku = ^ -— 7 ^ or? 

4 V - (1 + ^2)2 


Note that the l/5 2 -pole is spurious: its coefficient 
vanishes. 


( 14 ) 
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in the Brillouin zone give rise to sixteen propa¬ 
gating fermions. Correspondingly loop integrals 
which involve naive fermion propagators have six¬ 
teen soft regions. One of those is the region 
where all components of the loop momentum r]i 
are small. The fifteen additional doubler con¬ 
tributions arise after transforming one or sev¬ 
eral components of the integration momentum as 
r]i —> 1 /t]i and then expanding the resulting inte¬ 
grand around small 77 .;. For purely fermionic inte¬ 
grals, the sixteen soft contributions are all equal. 
In integrals with both boson and fermion prop¬ 
agators, the boson propagator is far off-shell for 
the doubler contributions and shrinks to a point 
upon expanding. 

The contribution of the hard part is obtained 
as usual, by expanding the loop integral in ex¬ 
ternal momenta and particle masses. The lattice 
tadpole integrals that occur in the case at hand 
are 


H({ai},n,m) = J] 


drji 


i =1 * 


(i 


■vf) 


1 


D 


,m+<5 ’ 


The regulator has to be on the fermion propaga¬ 
tor in order to regulate both the singularities at 
77 ^ = 0 and 77 ^ = 00 . The integration-by-parts 
and partial fractioning relations for this class of 
integrals are 


0 


1 + 2 (a* — 1) + 2 a i (a i — 1) n n 


T 2 cl 1 (CLi 1 ) (2 CLi 


1) (m + 5) m > H 


0 

0 


1 - ^2 a i (1 - a,i)m \ H , (15) 

i— 1 ' 

d 'I 

1 - ^(1 - a,i)n \ H . 


^ i= 1 J 

These relations are sufficient to reduce all inte¬ 
grals H to sixteen convergent master integrals. 


6. Summary and conclusion 

We have applied methods developed for the 
evaluation of continuum loop integrals to calcu¬ 
lations in lattice perturbation theory. The tech¬ 
nique of asymptotic expansions is used to expand 


lattice loop integrals around the continuum limit. 
After performing the expansion, all nontrivial de¬ 
pendence on momenta and masses is encoded in 
continuum loop integrals. The only genuine lat¬ 
tice integrals left are massless tadpole integrals. 
With the help of integration-by-part relations we 
then reduce all tadpole integrals to a small num¬ 
ber of master integrals. Except for the numerical 
evaluation of the master integrals, the entire cal¬ 
culation is performed analytically. 

The technique we have presented does not de¬ 
pend on a specific form of the lattice action and 
after illustrating it with a simple bosonic inte¬ 
gral, we have discussed how to apply it to the 
case of staggered fermions. Since the techniques 
we have been using were developed for multi-loop 
integrals, we hope that this method will also be 
useful beyond one loop. 

This research was supported by the DOE un¬ 
der grants number DE-AC03-76SF00515 and DE- 
FG02-97ER41022. 
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